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I.  ^n  anticipation  of  significant  access  to  the  NRL  CRAY  XMP-12 
supercomputer  we  devoted  major  attention  this  past  year  to  converting  our 
programs  to  that  computer,  which  runs  under  a  COS  operating  system.  We 
converted  the  following  programs  to  the  NRL  CRAY  XMP-12:  MOLASYS,  MRD-CI, 
CASSCF,  INTER-MOLASYS,  CRYSTAL  and  POLY-CRYST.  We  also  meshed  in  With  the 
MRD-CI  program  our  own  ab-initio  MOLASYS  MQDPOT/VfcDDO/MERGE  program  which 
contains  a  number  of  desirable  computational  strategies  for  ab-initio 
calculations  on  large  molecules  as  well  as  procedures  fer"Strictly 
orthogonal  localization  of  molecular  orbitals. 

II.  a  preliminary  to  carrying  out  MRD-CI  calculations  on  the  attack  of 
oxetane  on  protonated  oxetane  in  the  propagation  step  of  cationic 
polymerization,  we  initiated  MRD-CI  calculations  on  opening  the  oxetane  ring 
and  the  protonated  oxetane  ring.  These  calculations  were  carried  out  using 

strictly  orthogonal  localized  orbitals  in  the  region  of  the  0 — Cr^T  bond 

being  opened  with  the  remainder  of  the  occupied  molecular  orbitals  being 
folded  into  an  effective  Cl  Hamiltonian.  We  monitored  the  dominant 
configurations  in  the  final  wave  functions  as  a  function  of  the  Rq_^  bond 

distance.  We  also  contjmied^our  ab-initio  MODPOT  and  electrostatic 
molecular  potential  contour  calculations  around  substituted  oxetanes.  Among 
those  we  caTtrOTated  were  BNMAMO  and  MNMAMO. 

1 1 1. "He  carried  out  MRD-CI  calculations  for  the  dissociation  of  the  >C-N02 

bond  in  nitrobenzene.  Since  it  is  not  unambiguous  how  best  to  localize 
molecular  orbitals  in  an  aromatic  system  we  carried  out  the  calculations 
based  on  local  canonical  orbitals  in  the  region  of  the  >C-N02  portion  of  the 

nitrobenzene.  MRD-CI  calculations  were  carried  out  for  the  lowest  roots  of 

/  1313131  3 

all  eight  electronic  states:  A,,  A,,  A7,  A?,  B, ,  B,,  B?  and  B?  as  a 


3A,,  1A0,  3A„,  3B,,  and  3B, 


ail  ClyiiL  c  iclu  u  M  i  v  ouauca.  n^t  n  aiiu  ua  a 

function  of  C-N02  distance.  By  monitoring  the  dominant  configurations  of 

the  various  molecular  orbital  fillings  as  a  function  of  >C  -  N02  distance  it 

is  possible  to  determine  accurately  how  many  times  a  particular  potential 
energy  surface  changes  dominant  configuration.  A  nunber  of  the  higher 
electronic  states  were  predissociati ve.  These  predissociative  sta'tes  have 
significant  implications  for  the  initiation  of  detonation. 

IV.  He  continued  testing  of  our  POLY-CRYST  program  for  ab-initio 
calculations  in  crystals  and  polymers. 

V.  We  also  implemented  further  additions  to  the  INTER-MOLASYS  program 
including  the  capability  of  calculating  the  dispersion  energy  contributions 
by  a  variation-perturbation  procedure.  We  have  also  calculated  further  ab- 
initio  atom  class-atom  class  potential  functions.  \ 
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CONCISE  SUMMARY 


In  anticipation  of  significant  access  to  the  NRL  CRAY  XMP-12  supercom¬ 
puter  we  devoted  major  attention  this  past  year  to  converting  our  programs 
to  that  computer,  which  runs  under  a  COS  operating  system.  We  converted  the 
following  programs  to  the  NRL  CRAY  XMP-12:  MOLASYS  [Gaussian  integrals 
(including  special  options  for  large  molecules),  SCF,  Properf'  :>];  MRD-CI 
(Gaussian  Integrals,  Adapt,  SCF,  Transformations,  Multireference  double 
excitation.  Cl);  CASSCF  (Gaussian  integrals,  SCF,  Properties, 

Transformation,  MCSCF);  INTER-MOLASYS  (Integrals,  SCF,  Intermolecular  SCF, 
Special  routines  for  ab-initio  energy  partitioned  contributions  for  ab- 
initio  atom  class-atom  class  potentials  as  input  for  CRYSTAL);  CRYSTAL 
(optimization  of  crystal  unit  cell  parameters);  POLY-CRYST  [Crystal  and 
Polymer  SCF  and  integrals  (for  repeating  units].  We  also  meshed  in  with  the 
MRD-CI  program  our  own  ab-initio  MOLASYS  MODPOT/VRDDO /MERGE  program  which 
contains  a  number  of  desirable  computational  strategies  for  ab-initio  cal¬ 
culations  on  large  molecules  as  well  as  procedures  for  strictly  orthogonal 
localization  of  molecular  orbitals. 


As  a  preliminary  to  carrying  out  MRD-CI  calculations  on  the  attack  of 
oxetane  on  protonated  oxetane  in  the  propagation  step  of  cationic 
polymerization,  we  initiated  MRD-CI  calculations  on  opening  the  oxetane  ring 
and  the  protonated  oxetane  ring.  These  calculations  were  carried  out  using 


U 

strictly  orthogonal  localized  orbitals  in  the  region  of  the  0--C^![jj  bond 

being  opened  with  the  remainder  of  the  occupied  molecular  orbitals  being 
folded  into  an  effective  Cl  Hamiltonian.  We  carried  out  calculations  with 
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our  ab-initio  M00P0T  basis  set  and  with  an  all-electron  basis  set  contracted 

to  3S2*3.  We  monitored  the  dominant  configurations  in  the  final  wave  func¬ 
tions  as  a  function  of  the  Rq_q  bond  distance.  We  are  now  initiating  MRD-CI 

calculations  on  the  attack  of  oxetane  on  protonated  oxetane,  with  a  forma¬ 
tion  of  an  0 — C  bond  between  the  two  rings  and  opening  of  the  O-C^j  bond  in 

the  protonated  oxetane  ring.  These  calculations  are  being  carried  out  based 
on  strictly  orthogonal  localized  orbitals  in  the  region  extending  over  the 
new  bond  being  formed  and  the  opening  of  the  C-0^  bond.  We  also  continued 

our  ab-initio  MODPOT  and  electrostatic  molecular  potential  contour  calcula¬ 
tions  around  substituted  oxetanes.  Among  those  we  calculated  were  BNMAMO 
and  MNMAMO. 

We  carried  out  MRD-CI  calculations  for  the  dissociation  of  the  >C-N02 

bond  in  nitrobenzene.  Since  it  is  not  unambiguous  how  best  to  localize 
molecular  orbitals  in  an  aromatic  system  we  carried  out  the  calculations 
based  on  local  canonical  orbitals  in  the  region  of  the  >C-NQ2  portion  of  the 

nitrobenzene.  MRD-CI  calculations  were  carried  out  for  the  lowest  roots  of 
all  eight  electronic  states:  *Aj,  ^Ap  *A2,  ^A2,  *Bp  ^Bp  *B2  and  ^B2  as  a 
function  of  >C-N02  distance.  By  monitoring  the  dominant  configurations  of 
the  various  molecular  orbital  fillings  as  a  function  of  >C-N02  distance  it 
is  possible  to  determine  accurately  how  many  times  a  particular  potential 
energy  surface  changes  dominant  configuration.  The  ground  *Aj  state  is  not 
predissociati ve .  However,  the  relative  contributions  of  the  various  deter¬ 
minants  change.  The  ^  state  is  predissociative.  Both  the  *A2  and  ^A2 
states  are  predissociative  and  both  change  dominant  configuration  at  least 
twice  along  their  potential  energy  surfaces.  The  state  is  predissocia¬ 
tive  and  changes  dominant  configuration  at  least  once  along  its  potential 
energy  surface.  The  ^Bj  state  is  very  slightly  predissociative,  changing 
dominant  configuration  at  least  once  along  its  potential  energy  surface. 

The  *B2  and  3B2  surfaces  are  predissociative,  each  changing  dominant  con¬ 
figuration  at  least  twice  along  their  potential  energy  surfaces.  These 
results  emphasize  the  necessity  of  carrying  out  such  MRD-CI  calculations  for 
the  understanding  of  dissociation  of  nitroexplosives  and  initiation  of 
detonation . 

We  continued  testing  of  our  POLY-CRYST  program  for  ab-initio  calcula¬ 
tions  in  crystals  and  polymers. 


We  also  implemented  further  additions  to  the  INTER-MOLASYS  program 
including  the  capability  of  calculating  the  dispersion  energy  contributions 
by  a  variation-perturbation  procedure.  Our  testing  indicated  that  this 
latter  contribution  is  better  represented  using  the  virtual  orbitals  of  the 
two  interacting  molecules  than  by  just  using  the  virtual  orbitals  of  a 
single  molecule.  We  have  also  calculated  further  ab-initio  atom  class-atom 
class  potential  functions. 
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I.  Conversion  of  Our  Programs  to  CRAY  Vector  Supercomputers 


In  anticipation  of  significant  access  to  the  NRL  CRAY  XMP-12  supercom¬ 
puter  we  devoted  major  attention  this  past  year  to  converting  our  programs 
to  that  computer,  which  runs  under  a  COS  operating  system.  The  other  CRAY 
computers  on  which  we  ran  up  until  this  year  have  all  had  various  CTSS 
operating  systems.  Our  programs  are  large  and  complicated  and  also  make  use 
of  the  special  mathematics  libraries  on  the  computers.  None  of  the  CRAY 
computers  we  use  are  completely  compatible  and  thus  we  have  to  adapt  the 
programs  separately  for  each  computer. 

We  began  converting  to  the  NRL  CRAY  XMP-12  in  March  1985  just  as  soon 
as  that  computer  was  available.  We  have  had  helpful  cooperation  from  the 
NRL  computer  center  and  also  interacted  closely  with  them  in  testing  exter¬ 
nal  data  transmission. 

We  had  expected  that  the  NRL  CRAY  XMP-12  would  be  our  major  production 
computer.  Thus  all  of  our  programs  would  be  under  one  operating  system 
(COS)  which  is  vital  ,  since  many  of  these  large  programs  have  to  interface 
with  our  other  large  programs.  For  some  of  the  programs,  assembly  language 
routines  have  been  implemented  to  make  the  handling  of  large  data  files 
(integrals,  transformed  integrals,  etc.)  much  more  tractable.  These  as¬ 
sembly  language  portions  were  written  in  CAL,  the  CRAY  assembly  language 
running  under  COS.  CAL  is  not  available  on  the  CRAY’S  which  run  under  CTSS. 

We  converted  the  following  programs  to  the  NRL  CRAY  XMP-12 

MOLASYS  -  Gaussian  integrals  (including  special  options  for 
large  molecules ) 

-  SCF 

-  Properties 

MRD-CI  -  Gaussian  Integrals 

-  Adapt 

-  SCF 

-  Transformations 

-  Multireference  double  excitation  -  Cl 

Earlier  in  this  past  year  before  the  NRL  CRAY  XMP  became  available  we 
had  converted  the  complete  MRD-CI  program  to  a  CRAY  1M  (running  under  CTSS). 
This  had  involved  a  great  deal  of  rewriting.  Then  when  the  NRL  CRAY  XMP 
became  available  we  made  the  MRD-CI  program  operational  on  the  NRL  CRAY  XMP. 


The  NRL  machine  has  sufficient  available  disc  space  to  make  feasible 
reasonably  large  MRD-CI  calculations  on  good  size  molecules.  We  have  been 
making  extensive  use  of  the  MRD-CI  program  both  for  decomposition  of  ener¬ 
getic  molecules  and  most  recently  for  the  propagation  step  for  cationic 
polymerization.  We  also  meshed  in  with  the  MRD-CI  program  our  own  ab-initio 
MOLASYS  MODPOT/VRDDO/MERGE  program  which  contains  a  number  of  desirable 
computational  strategies  for  ab-initio  calculations  on  large  molecules  as 
well  as  procedures  for  strictly  orthogonal  localization  of  molecular 
orbitals. 

CASSCF  -  Gaussian  integrals 

-  SCF 

-  Properties 

-  Transformation 

-  MCSCF 

We  have  also  vectorized  the  computer-time-consuming  steps  of 
this  latter  series  of  programs.  All  of  the  I/O  routines  had  to  be 
rewritten  to  run  under  COS. 

INTER-MOLASYS  -  Integrals 

-  SCF 

-  Intermolecular  SCF 

-  Special  routines  for  ab-initio  energy  partitioned 
contributions  for  atom  class-atom  class  potentials  as 
input  for  CRYSTAL 

CRYSTAL  -  Optimization  of  crystal  unit  cell  parameters 

POLY-CRYST  -  Crystals  and  polymer  SCF  and  integrals  (for  repeating 
units ) 

Our  ONR  technical  contract  monitor,  Dr.  Richard  Miller,  had  even  hoped 
that  it  would  be  possible  to  arrange  high  speed  55  kilobit  data  transmission 
to  and  from  the  NRL  CRAY  to  our  lab  at  the  Johns  Hopkins  University.  Thus 
we  researched  intensively  the  network  and  equipment  facilities  necessary  for 
such  transmission.  We  presented  to  our  ONR  technical  contract  monitor  in 
July  1985  the  complete  system  hardware  and  software  requirements,  procedures 
and  persons  to  contact  in  DOD  Telecommunications.  However,  in  view  of  the 
fact  that  our  ONR  technical  contract  monitor  and  we  learned  only  at  that 
same  time  of  our  very  limited  computer  time  allottment  on  the  NRL  CRAY  XMP- 
12,  the  high  speed  data  transmission  project  is  on  hold  for  the  time  being. 

However,  in  mid  July  1985  unexpectedly  we  were  informed  that  only  a  few 
hours  of  computer  time  would  be  allotted  to  us  on  the  NRL  CRAY  XMP-12. 
Moreover,  the  NRL  CRAY  computer  time  which  was  to  have  been  available  at  no- 
cost  up  until  October  1,  1985  began  being  charged  (partial  rates)  as  of 
August  1,  1985. 

We  have  been  very  fortunate  to  get  computer  time  on  other  CRAY 
computers: 

CRAY  1M  -  CTSS 

CRAY  IS  -  COS  [but  without  the  special  NAG  (Nottingham 
Algorithm  Group)  library] 
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Thus,  we  have  also  converted  some  of  our  programs  to  these  CRAY 
computers.  Our  programs  are  large  and  complicated  and  require  considerable 
effort  each  time  we  move  to  a  different  CRAY  since  no  two  of  the  machines 
are  completely  compatible  systemwise. 

This  year  we  have  been  carrying  out  production  runs  on  various  CRAY’S 
with  a  number  of  our  programs. 
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II.  Calculations  for  Cationic  Polymerization  of  Energetic  Monomers 

The  Navy  has  an  interest  in  energetic  polymers  made  by  cationic 
polymerization  from  cyclic  ethers  substituted  with  exotic  energetic  pendant 
groups . 

Previous  experimental  results  of  Aoki  et  al .  on  cationic  polymerization 
of  cyclic  ethers  had  indicated  that  the  ease  of  polymerization  was  related 
to  the  basicity  of  the  cyclic  ether  and  to  its  ring  strain.  By  generation 
of  electrostatic  molecular  potential  contour  (EMPC)  maps  from  ab-initio 
M00P0T/VRDD0/MER6E  quantum  chemical  calculations  we  were  able  to  predict 
correctly  the  order  of  the  basicity  of  oxetanes  substituted  with  exotic 
energetic  groups  even  prior  to  their  actual  systhesis.  Our  EMPC  maps  also 
predicted  correctly  the  order  of  the  polymerizabil ity  of  these  substituted 
oxetanes  even  prior  to  the  polymerization  experiments.  Also  valuable  was 
our  prediction  from  the  EMPC  maps  that  3,3-dinitrooxetane  should  not 
polymerize  under  ordinary  conditions.  Experimentally  this  was  shown  to  be 
so.  Thus,  we  now  even  have  a  "fingerprint"  EMPC  map  which  can  be  compared 
to  EMPC  maps  calculated  even  for  as  yet  unsynthesized  substituted  oxetanes 
to  predict  whether  or  not  that  particular  oxetane  will  polymerize. 

Our  calculations  of  the  ring  strain  for  variously  substituted  oxetanes 
showed  little  variation  of  ring  strain  with  different  substitutents .  Thus 
the  change  in  the  basicity  of  the  oxetanes  with  substitution  was  the  over¬ 
riding  factor  in  determining  their  polymerizabil ities . 

The  propagation  step  had  been  suggested  earlier  by  experimentalists  to 
be  influenced  by  the  electrophil icity  of  the  resulting  ring-opened 

\/Y 

/c\ 

"protonated  carbocation" .  For  a  substituted  oxetane  >C  C<,  the  ring- 

X 

i  + 

opened  "protonated  carbocation"  would  be  H0-CHo-C-CHo.  There  are  several 

*  i  £■ 

Y 

calculated  quantum  chemical  indices  which  relate  to  electrophil icity:  the 
charge  on  the  carbon  atom  of  the  terminal  -CH^,  the  EMPC  map  around  that 

terminal  -CH^  and  the  energies  of  the  unoccupied  molecular  orbitals  of  the 

"protonated  carbocation."  We  had  calculated  all  of  these  indices.  There 
was  little  difference  in  the  charges  on  the  carbon  atoms  on  the  terminal 

-CH^  group  or  in  the  EMPC  maps  around  that  group  as  a  function  of 

substitution.  However,  there  was  a  significant  difference  in  the  energies 
of  the  unoccupied  molecular  orbitals  as  a  function  of  substitution.  On  this 
latter  criterion  the  ring-opened  "protonated  carbocation"  of  3,3- 
dinitrooxetane  was  the  most  electrophilic  species.  However,  3,3- 
dinitrooxetane  will  not  polymerize  under  ordinary  conditions  because  its 
basicity  is  too  low.  Therefore,  it  is  the  initiation  step  rather  than  the 
propagation  step  which  is  overriding  in  determining  the  polymeri zabi 1 ity  of 
a  monomer . 


The  next  aspect  of  the  mechanism  of  cationic  polymerization  is  the 
relative  polymerization  ratio  in  copolymeri zation . 

We  have  been  exploring  what  quantities  to  calculate  to  predict  relative 
polymerization  ratios.  This  turns  out  to  be  considerably  more  complicated 
than  initiation.  We  can  predict  correctly  in  a  mixture  which  monomer  will 
initiate  first  and  more  quickly.  However,  propagation  is  dependent  on  a 
number  of  different  and  sometimes  conflicting  factors:  the  electrophil icity 
of  the  resulting  ring-opened  "protonated  carbocation"  (which  substituents 
influence  in  just  the  reverse  manner  to  how  they  influence  initiation),  the 
accessibility  of  the  negative  EMPC  map  around  the  ring  oxygen  (the  EMPC  map 
is  related  to  the  basicity  of  the  oxygen)  and  steric  hindrance  to  attack  on 
the  ring  oxygen.  There  have  been  a  few  experimental  studies  on  relative 
polymerization  ratios  in  this  overall  Navy  f-^ject,  but  to  date  these 
partners  have  been  picked  for  their  relevance  to  the  actual  polymers  being 
synthesized.  However,  these  monomers  have  not  been  picked  to  make  a 
coherent  series  for  comparison  to  theory.  We  have  had  continuing  fruitful 
detailed  discussions  with  Gerald  Manzer,  who  is  synthesizing  a  number  of 
monomers  and  polymers.  A  few  well-defined  experiments  could  help  disen¬ 
tangle  some  of  the  conflicting  and  complicating  factors  in  the  propagation 
step.  Manzer  indicated  to  us  in  June  and  again  in  recent  conversations  that 
he  will  be  carrying  out  some  fundamental  studies  of  cationic  copolymeriza¬ 
tion  with  standards  such  as  oxetane  to  have  a  set  of  coherent  results  for 
meaningful  correlations  of  copolymerization  ratios. 

In  a  copolymerization  experiment  monomers  1  and  2  must  first  undergo 
initiation.  Hence,  for  a  preliminary  attempt  to  compare  our  calculated 
theoretical  indices  with  experimentally  measured  copolymerization  reactivity 
ratios,  we  considered  our  EMPC  maps,  which  are  indicative  of  basicity  and 
hence  of  propensity  of  a  monomer  to  initiate  (or  preference  of  a  monomer  to 
be  attacked  by  a  propagating  species). 

Formation  of  random  copolymers  are  determined  by  the  selectivity  of  the 
propagation  reaction,  i.e.  the  relative  reactivities  of  the  comonomers. 
However,  it  has  been  shown  experimentally  for  copolymerization  of  more 
conventional  species  that  relative  monomer  reactivities  on  the  copolymeriza¬ 
tion  may  differ  strongly  from  the  apparent  monomer  reactivities  in 
homopolymerization.  A  species  which  is  less  reactive  in  homopolymerization 
may  be  the  preferentially  consumed  comonomer  in  the  cationic 
copolymerization.  This  complicates  greatly  the  prediction  of  relative 
polymerization  ratios  for  copolymerization. 

We  are  also  considering  some  additional  theoretical  approaches  to 
understanding  the  propagation  step. 

Recently,  Gerry  Manzer,  Dr.  Lillya  and  Dr.  Bob  Earl  have  suggested  to 
us  that  rather  than  a  ring-opened  "protonated  carbocation",  more  direct 
attacks  might  take  place. 

The  propagation  step  might  be  more  complicated  than  ring  opening  of  the 
"protonated  carbocation"  and  subsequent  attack  of  the  ring  opened 
"protonated  carbocation"  on  the  second  molecule  of  monomer.  Rather,  they 
hypothesized  that  there  might  be  an  attack  of  the  "protonated  carbocation" 
on  the  second  molecule  of  monomer  with  a  concommitant  ring-opening. 
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the  mechanism  of  the  propagation  step  in  cationic  polymerization.  Both  ab- 
initio  M0DP0T  and  ab-initio  all-electron  valence  double  zeta  MRD-CI 
calculations  were  carried  out  to  investigate  basis  set  dependence.  The 
results  indicated  that  the  potential  energy  surfaces  resembled  each  other 
whether  only  localized  orbitals  in  the  region  of  the  bond  being  broken  or 
whether  the  full  set  of  canonical  molecular  orbitals  were  used  in  the  MRD- 
CI.  So  far  the  MRD-CI  potential  energy  surfaces  have  been  investigated  with 
ab-initio  M0DP0T  or  an  ab-initio  all-electron  valence  double  zeta  basis 
sets.  Further  basis  sets  will  be  investigated  but  these  initial  results  are 
very  encouraging.  The  preliminary  all -electron  MRD-CI  results  indicate  that 
it  takes  less  energy  to  open  the  protonated  oxetane  ring  than  the  ring  of 
neutral  oxetane.  We  are  now  carrying  out  preliminary  ab-initio 
M0DP0T/VRDD0 /MERGE  MRD-CI  calculations  for  opening  the  protonated  oxetane 
ring  in  the  presence  of  another  attacking  oxetane  molecule. 


This  investigation  uses  our  strategy  of  explicit  MRD-CI  calculations  on 
the  localized  orbitals  in  the  region  of  the  bond  with  the  remainder  of  the 
occupied  molecular  orbitals  being  folded  into  an  effective  Cl  Hamiltonian. 

This  study  of  opening  the  protonated  oxetane  ring  while  a  second 
oxetane  molecule  is  attacking  is  an  excellent  illustration  of  the  power  and 
versatility  of  the  various  programs  which  we  have  meshed  together  and  gotten 
operational  on  the  CRAY'S.  Without  using  our  ab-initio  M0DP0T/VRDD0 /MERGE 
procedure  the  integral  times  would  become  very  large  for  higher  X  and  Y 
substituents  and  for  the  number  of  geometries  that  would  have  to  be 
investigaged .  Without  carrying  out  the  MRD-CI  calculations  based  on  local¬ 
ized  orbitals  in  the  region  of  the  bond  being  opened,  the  MRD-CI 
calculations  would  be  intractable  in  this  generation  of  CRAY  1  or  CRAY  XMP 
computers.  However,  the  localized  orbital  MRD-CI  approach  looks  extremely 
promising  for  investigating  the  propagation  step  of  the  substituted 
oxetanes.  With  this  technique  we  should  be  able  to  gain  insight  into  such 

questions  as  will  the  propagation  step  in  copolymerization  (where  X  and  X' 

are  different  and  Y  and  Y'  are  different)  prefer  to  go 
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This  same  series  of  computer  programs  is  also  being  and  will  continue 
to  be  used  for  bond  breaking  problems  on  dissociation  of  energetic  molecules 
and  as  building  blocks  for  a  sequence  of  programs  to  calculate  input  data 
for  fractoemission  calculations  and  for  calculations  on  fractoemission 
processes . 


We  had  originally  intended  to  carry  out  the  calculations  for  this 
oxetane  propagation  step  investigation  on  the  NRL  CRAY  XMP-12  as  soon  as  we 
got  all  of  the  necessary  program  pieces  converted  and  meshed  together.  The 
very  initial  studies  were  carried  out  on  the  NRL  CRAY  XMP-12.  However  in 
view  of  the  severely  limited  amount  of  computer  time  allotted  to  us  on  the 
NRL  CRAY,  we  had  to  spend  time  again  to  convert  the  entire  meshed  program 
sequence  to  another  CRAY  where  a  more  realistic  amount  of  computer  time  is 
available.  This  conversion  was  accomplished  and  the  calculations  are  well 
underway. 


AB-INITIO  MULTIREFERENCE  DOUBLE  EXCITATION  CONFIGURATION  INTERACTION 

MRO-CI 

1.  DO  A  MODEST  Cl  CALCULATION  (SEVERAL  THOUSAND  CONFIGURATIONS). 

2.  FROM  THESE  PICK  THE  MOST  IMPORTANT  (USUALLY  -  20). 

3.  ALLOW  ALL  SINGLE  AND  DOUBLE  EXCITATIONS  FROM  ALL  REFERENCE 
DETERMINANTS  (GENERATES  -  750,000  CONFIGURATIONS). 

4.  ESTIMATE  THE  ENERGY  CONTRIBUTION  OF  EACH  OF  THE  -  750,000 
CONFIGURATIONS  BY  A  PERTURBATIVE  PROCEDURE. 

5.  SET  A  THRESHOLD.  DO  Cl  EXPICITLY  FOR  CHOSEN  CONFIGURATIONS. 

6.  ADD  BACK  IN  THE  ENERGIES  OF  OTHER  CONGIGURATIONS  BY  PERTURBATIVE 
PROCEDURE  AND  EXTRAPOLATE. 

7.  ADD  DAVIDSON  TYPE  CORRECTION  FOR  SIZE  CONSISTENCY. 


E(FULL  Cl)  “  E(EXT)  +  (1  ~  l  CP  ^E(EXT)  "  E ( RE F ) ^ } 


The  following  tables  present  examples  of  the  types  of  information  which 
result  from  our  MRD-CI  calculations  on  opening  the  oxetane  ring. 


Hn  ,H 

Hs  C^  H 

H^CC 

H  R^O^  SH 

The  MRD-CI  calculations  were  run  from  R  -  2.8  a.u.  [equilibrium 
geometry  -  the  undistorted  oxetane  (or  protonated  oxetane)  ring]  to  4.9  a.u. 
[the  opened  oxetane  (or  protonated  oxetane)  ring]. 


•  Jfc  *  4  »  * 
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Tables  I-IV  present  the  number  of  active  orbitals  in  the  MRD-CI,  which 
specific  localized  molecular  orbitals  were  included  in  the  MRD-CI,  the  ES(,p, 

Egxi*  EfuU  and  the  coefficients  (c1)  of  the  dominant  configurations  in  the 

final  MRD-CI  wave  function.  For  these  initial  test  calculations,  both  our 
ab-initio  MODPOT  basis  was  used  (valence  electrons  only)  as  well  as  a  more 

conventional  Dunning  3S3^  contraction  of  a  Huzinaga  9S5^  atomic  basis  set. 


Table 

I 

Oxetane 

(R  ■  2.8  a .u . ) 

MRD-CI 

Table 

II 

Oxetane 

(R  -  4.9  a.u.) 

MRD-CI 

Table 

III 

Protonated  Oxetane 

(R  *•  2.8  a.u.) 

MRD-CI 

Table 

IV 

Protonated  Oxetane 

(R  -  4.9  a.u.) 

MRD-CI 

Tables  V-VI 1 1  present  the  energies  (EgQp  *  Eqj.  and  Efuv|)  and  the 

dominant  configurations  in  the  final  MRD-CI  wave  functions  of  oxetane  and 
protonated  oxetane  as  a  function  of  R. 


OXETANE  (equilibrium  geometry  R  -  2.8  a.u.)  MRD-CI  Energies  (a. 
MODPOT  Basis  Dunning  (3s2p)  basis 
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This  past  year  we  also  carried  out  some  additional  quantum  chemical 
calculations  and  generated  EMPC  maps  for  substituted  oxetanes.  The  oxetanes 
we  investigated  were: 

BMNAMQ  3,3-bis(methylnitrami nomethyl  )oxetane 

for  both  the  experimental  crystal  structure,  which  had 
a  puckered  oxetane  ring,  and  for  a  planar  oxetane  ring, 
which  was  the  lower  in  calculated  total  energy. 

MMNAMO  3-methyl -3(methylni traminomethyl )oxetane 

for  both  a  puckered  oxetane  ring  (as  in  the  crystal 
structure  of  MNBMNAMO)  and  for  a  planar  oxetane  ring, 
which  again  was  lower  in  total  calculated  energy. 

We  will  still  be  carrying  out  more  calculation  for  optimization  of 
geometries  of  these  molecules  since  these  compounds  have  multiple  maxima  and 
minima . 

These  oxetanes  with  the  floppy  substituents  have  multiple  maxima  and 
minima.  Thus  merely  optimizing  geometries  with  any  customary  derivative 
method  can  only  lead  to  a  local  minimum  dependent  on  the  starting  geometry. 
We  plan  to  give  attention  to  developing  a  program  strategy  which  will  com¬ 
bine  our  ab-initio  MODPOT/VRDDO/MERGE  technique  with  derivatives.  Our 
recent  MRD-CI  calculations  on  oxetane  and  protonated  oxetane  with  both  ab- 
initio  MODPOT  and  ab-initio  all  electron  double  zeta  atomic  basis  sets 
indicate  the  reliability  of  the  ab-initio  MODPOT  atomic  basis  sets  which  can 
be  used  for  the  vast  majority  of  geometry  optimizations  of  oxetanes. 
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III.  MRD-CI  Calculations  on  Decomposition  Pathways  of  Energetic 

Nitrocompounds 

A.  Nitrobenzene 

The  dissociation  of  >C  -  NO2  and  >N  -  NO^  bonds  of  nitroexplosi ves  are 

one  of  the  key  primary  steps  in  initiation  of  explosives.  It  had  been 
suggested  recently  that  in  nitroexplosi ves  where  the  potential  energy  sur¬ 
faces  for  ground  or  electronically  excited  states  indicated  that  one  or  more 
of  the  dissociative  pathways  could  undergo  predissociation,  these  predis- 
sociative  pathways  could  lead  to  products  with  sizable  amounts  of  kinetic 
energy,  thus  initiating  processes  leading  to  detonation  in  the  crystal.  The 
molecular  decomposition  pathways  of  the  ni troexplosi ves  and  the  initiation 
and  subsequent  steps  in  detonation  would  then  be  completely  intertwined.  To 
study  molecular  decomposition  pathways  it  is  necessary  to  use  ab-initio  MRD- 
CI  (multireference  determinant-conf iguration  interaction)  or  CASSCF 
(complete  active  space  multiconfiguration  SCF)  calculations.  Our  previous 
MRD-CI  and  CASSCF  results  on  the  >C  -  N02  decomposition  of  nitromethane  and 

our  MRD-CI  calculations  on  the  >N  -  NO^  decomposition  pathway  of  RDX  had 

indicated  that  there  were  multiconfiguration  electronic  wave  functions  even 
in  the  ground  electronic  state  at  equilibrium  geometry  and  additional  con¬ 
figurations  along  the  dissociation  pathway.  Because  the  nitroexplosi ves  are 
large  molecules,  it  is  still  computationally  rather  intractable  computer 
space  and  time  wise  to  carry  out  Cl  calculations  for  the  entire  molecule. 

By  localizing  and  focusing  on  the  >C  -  NO2  or  >N  -  N02  decomposition  path¬ 
way,  much  important  insight  can  be  gained  about  how  many  and  what  kind  of 
potential  energy  surfaces  are  involved,  and  to  what  types  of  product  species 
they  lead.  This  MRD-CI  method,  which  also  includes  techniques  to  account 
for  linked  double  excitations,  is  also  applicable  to  intermolecular 
interactions.  The  most  rational  strategy  on  this  generation  of  CRAY-1  or 
CRAY-XMP  supercomputers  is  to  carry  out  the  MRD-CI  or  CASSCF  calculations  on 
localized  or  local  molecular  orbitals  in  the  region  of  the  bond  being  broken 
(either  strictly  orthogonal  localized  orbitals  or  canonical  orbitals  in  that 
region  of  space)  with  the  remainder  of  the  occupied  molecular  orbitals  being 
folded  into  an  effective  Cl  Hamiltonian. 

For  our  initial  studies  on  the  >C  -  N02  decomposition  pathway  of 

nitrobenzene  we  carried  out  ab-initio  all-electron  MRD-CI  calculations  with 

a  10s 6P  -*•  3S2P  basis  set  on  the  C,  N,  0,  atoms  and  a  3s  -*■  Is  basis  on  the  H 
atoms.  There  were  a  total  of  50  molecular  orbitals  (32  occupied  and  18 
virtual).  While  we  have  carried  out  ab-initio  strictly  orthogonal  localized 
orbital  MRD-CI  calculations  for  the  bond  dissociation  of  aliphatic  nitrocom¬ 
pounds  (such  as  RDX),  it  is  not  completely  unambigous  how  to  best  to 
localize  molecular  orbitals  for  an  aromatic  system.  Thus  we  carried  out  the 
MRD-CI  calculations  on  nitrobenzene  using  canonical  molecular  orbitals 
taking  explicit  account  in  the  MRD-CI  of  the  occupied  molecular  orbitals  in 

this  region  of  the  ^>C  -  NC^  bond  being  dissociated  (15  occupied  molecular 

orbitals)  and  all  of  the  virtual  orbitals.  The  other  17  of  the  occupied 
molecular  orbitals  were  folded  into  an  effective  Cl  Hamiltonian. 


Table  TX  presents  for  the  ^2»  ^2,  ^  and  ^ 

states  of  nitrobenzene  the  number  of  reference  configurations,  the  total 
number  of  configurations  generated,  the  number  of  configurations  selected 
(by  the  perturbative  procedure  with  a  threshold  of  30  pH). 
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Table  IX 


NITROBENZENE  MRD-CI  NUMBER  OF  CONFIGURATIONS  (RC_N»  2-8  A.U.) 


STATE 

NO.  OF  REFERENCE 
CONFIGURATIONS 

NO.  OF  CONFIGURATIONS 
GENERATED 

NO  OF  CONFIGURATIONS 
SELECTED 

\  (GS) 

8 

40662 

2651 

17 

576557 

2094 

lft2 

17 

98609 

2060 

3a2 

17 

98604 

2094 

1bi 

16 

90255 

2040 

\ 

8 

45748 

2010 

‘B2 

20 

111563 

2038 

16 

87816 

2106 

THRESHOLD  30pH 


As  with  the  dissociation  of  other  >C  -  N02  or  >N  -  N02  compounds,  there 

is  a  crossing  of  molecular  orbitals  along  the  dissociation  pathway  (Figure 
1).  The  surfaces  for  E^j,  E^y  and  are  all  essentially  parallel 

(Figure  2) . 

Our  MRD-CI  results  on  the  3A^  3A2,  *B2  and  3B2 

states  of  nitrobenzene  show  a  wealth  of  structure  in  the  potential  energy 
surfaces  for  the  various  electronic  states  of  nitrobenzene  as  a  function  of 
>C  -  N02  distance.  (Figure  3) 
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It  is  obvious  from  examination  of  the  potential  energy  surfaces  in  H 

Figure  3  that  there  are  a  number  of  potential  energy  surfaces  which  exhibit  1 

surface  crossing  behavior  and  some  which  even  exhibit  predissociative  I 

behavior.  To  verify  these  apparent  surface  crossings,  the  c*  for  the  \ 

dominant  configurations  in  the  final  MRD-CI  calculation  for  each  different  j 

surface  at  each  are  presented  in  Tables  X-XVII.  Examination  of  Table  t 

X  (the  ground  electronic  state  *Aj)  shows  that  even  at  equilibrium  geometry  I 

there  is  still  the  contribution  of  the  excitation  fj 

(if) *  -  (if*) *C(la* )*  -  (4b,)*] 

v 

in  the  NO,  group  (as  in  nitromethane) .  Also,  the  results  in  Table  X  verify  r 

that  a  number  of  other  configurations  become  important  along  the  dissocia-  I 

tion  pathway  even  of  the  ground  electronic  state.  For  some  of  the  higher  r 

states  which  exhibit  surface  crossings  it  can  be  seen  that  at  points  cor-  > 

responding  to  surface  crossings  there  is  a  complete  change  in  the  dominant  \ 

configurations  in  the  wave  functions. 

V* 

By  monitoring  the  dominant  configurations  of  the  various  molecular  M 

orbital  fillings  as  a  function  of  >C  -  N02  distance  it  is  possible  to  deter-  * 

mine  accurately  how  many  times  a  particular  potential  energy  surface  changes  / 

T 

dominant  configuration.  The  ground  A^  state  is  not  predissociative.  K 

However,  the  relative  contributions  of  the  various  determinants  change  as  a  |j 

function  of  >C  -  N02  distance  (Table  X).  The  3Aj  state  is  predissociative  >J 

13  !J 

(Table  XI).  Both  the  A2  and  A2  states  are  predissociative  and  both  change  /j 

dominant  configuration  at  least  twice  along  their  potential  energy  surfaces.  \ 

(Tables  XII  and  XIII)  | 

The  *B^  state  is  predissociative  and  changes  dominant  configuration  at 

least  once  along  its  potential  energy  surface  (Table  XIV).  > 

S 

The  ^Bj  state  is  very  slightly  predissociative,  changing  dominant  | 

configuration  at  least  once  along  its  potential  energy  surface  (Table  XV).  !- 

.* 

13 

The  B,  and  B?  surfaces  are  predissociative,  each  changing  dominant  i 

configuration  at  least  twice  along  their  potential  energy  surfaces  (Tables  k 

XVI  and  XVII).  \ 
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Thus  nitrobenzene  and  other  nitroaromatic  explosives  will  have  a  multi¬ 
tude  of  possible  predissociati ve  pathways. 

Hence,  in  summary,  it  is  apparent  that  detailed  MRD-CI  calculations 
such  as  we  report  here  will  be  necessary  to  describe  the  properties  and 
dissociation  of  both  ground  and  excited  electronic  states  of 
nitroexplosives .  This  will  be  especially  crucial  in  understanding  initia¬ 
tion  of  detonation  as  well  as  subsequent  steps. 

We  are  continuing  these  studies  to  identify  the  structures  of  the 
higher  lying  potential  energy  surfaces.  In  order  to  calculate  the  higher 
roots,  the  program  dimensions  will  have  to  be  enlarged  to  handle  a)  more 
reference  configurations,  b)  more  active  electrons  and  c)  to  pull  even  more 
higher  eigenvalues  from  the  large  Cl  matrix  than  the  program  now  does.  The 
increase  in  these  dimensions  will  be  important  for  research  on  any  larger 
energetic  compounds. 
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IV.  POLY-CRYST 

The  POLY-CRYST  program  (for  ab-initio  calculations  on  crystals  and 
polymers)  has  been  converted  to  the  CRAY-1M  and  to  the  CRAY-XMP.  We  in¬ 
itiated  preliminary  timing  test  calculations  on  a  system  of  three  unit  cells 
of  nitromethane  (four  molecules  per  unit  cell).  As  with  the  quantun  chemi¬ 
cal  calculations  of  many  nitrocompounds,  the  crystal  orbital  SCF  calculation 
of  nitromethane  is  difficult  to  converge.  Using  the  special  extrapolation 
and  damping  techniques  we  had  developed  over  the  years,  we  converged  this 
calculation.  However,  even  on  a  CRAY-XMP,  calculation  of 
the  three  unit  cell  system  on  nitromethane  is  not  trivial  timewise,  even 
using  a  minimal  atomic  basis  set,  ab-initio  effective  core  model  potentials 
for  the  inner  shells  and  our  charge  conserving  integral  prescreening 
evaluation.  The  calculation  for  one  unit  cell  of  nitromethane  (integrals 
and  SCF)  took  120  seconds.  For  the  three  cell  system  the  integrals  took  400 

seconds,  and  the  complex  SCF  for  20  iterations  but  only  six  £  points  (since 
the  three  cells  were  linear)  took  1100  seconds.  We  will  systematically  be 
increasing  the  number  of  unit  cells.  At  the  least,  27  unit  cells  should  be 
run  for  this  problem.  For  a  three-dimensional  system  one  may  have  to  use 

50-100  £  points.  We  will  be  investigating  this  as  well  as  whether  long 
range  effects  have  to  be  included  for  molecular  crystals  and,  if  so,  the 
best  way  to  include  them. 

Because  of  the  size  of  the  ni troexplosi ve  molecules  and  the  large  number 
of  molecules  in  a  unit  reference  cell  (as  high  as  8  for  RDX),  the  full  power 
of  our  POLY-CRYST  program  for  nitroexplosi ves  will  only  reach  fruition  on 
the  next  generation  of  CRAY-2  type  supercomputers.  (CRAY-2  computers,  which 
have  a  CPU  core  memory  of  256  million  real  64-bit  words,  have  already  been 
delivered.  This  compares  to  the  only  two  million  to  eight  million  words  of 
CPU  memory  on  the  current  CRAY-XMP 's.)  We  are  poised  and  ready  for  the  next 
generation  of  supercomputers. 
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V.  Ab-initio  Atom  Class  -  Atom  Class  Potential  Functions  and  Prediction  of 
Crystal  Densities 

In  addition  to  being  useful  in  predicting  optimum  crystal  unit  cell 
parameters  for  energetic  molecules,  these  ab-initio  atom  class  -  atom  class 
potential  functions  can  be  used  to  gain  insight  into  the  bonding  and  ar¬ 
rangement  of  such  molecules  in  clusters  and  the  interactions  of  such 
molecules  both  in  gas  and  condensed  phases.  Such  potentials  can  provide 
valuable  input  for  molecular  dynamics  and  hydrodynamics  modeling. 

We  converted  to  the  CRAY-XMP  our  programs  necessary  for  these  studies: 

INTER-MOLASYS  for  ab-initio  intermolecular  SCF  calculations 
with  energy  partitioning  to  obtain  atom 
class-atom  class  ab-initio  potential  energy 
functions  and  correction  for  basis  set  super¬ 
position  error  ( BSSE) . 


CRYSTAL-JHU  for  optimizing  unit  cell  parameters. 

We  have  derived  and  implemented  procedures  to  calculate  ab-initio  atom 
class-atom  class  potentials  from  energy  partitioned  ab-initio  SCF  calcula¬ 
tions  between  molecules  and  to  use  these  potentials  plus  dispersion  energy 
to  optimize  crystal  unit  cell  parameters  (from  which  densities  are 
calculated) . 


There  are  a  number  of  terms  in  our  expression  for  the  intermolecular 
energy  in  a  crystal.  The  main  components  of  the  two-body  interaction  energy 
are 


E  .  E(1>  ♦  E<1}  ♦  E(2)  ♦  E(2) 
t  tEL  fcEX  tIND  tDISP 


The  first  three  contributions  to  equation  I  can  be  obtained  from  the 


decomposition  of  the  variational  E^p 


interaction  energy  corrected  for  basis 


set  superposition  error  (BSSE). 


We  had  also  developed  procedures  to  calculate  the  E^  Mjp  (the  electros¬ 
tatic  contribution  from  the  atomic  multipole  moments)  and  the 
E inq  le  (the  long  range  induction)  explicitly  each  time  for  each  arrangement 

of  molecules.  The  short  range  terms  E^  pEN,  Eq_t,  Ej^  ct  are  fitted  to 

functional  forms.  The  dispersion  terms  are  estimated  semi-theoretical ly. 

We  initiated  considering  additional  approaches  to  calculating  these  poten¬ 
tials  (especially  the  dispersion  contributions  which  contribute  a 
significant  fraction  of  the  intermolecular  interaction  energy  even  for 
nitroexplosives  which  also  have  strong  charge  redistribution  and  hence 
strong  electrostatic  and  inductive  contributions  to  the  intermolecular 
interaction  energy). 
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